January 22, 2013 

FRACTIONAL ABSOLUTE MOMENTS OF HEAVY TAILED DISTRIBUTIONS 



MUNEYA MATSUI AND ZBYNEK PAWLAS 

Abstract. Several convenient methods for calculation of fractional absolute moments are given 
with application to heavy tailed distributions. We use techniques of fractional differentiation to 
obtain formulae for E[\X — /i| 7 ] with 1 < 7 < 2 and /i £ R. The main focus is on heavy tailed 
distributions, several examples are given with analytical expressions of fractional absolute moments. 
As applications, we calculate the fractional moment errors for both prediction and parameter esti- 
mation problems. 



1. Fractional moments 

In many diverse areas an important task is to predict future random elements based on previous 
available observations. For linear time series models, e.g. ARMA, ARIMA or FARIMA models, 
linear predictors given by linear combinations of past observations are often considered. More 
sophisticated statistical methods are used for prediction in non-linear models, such as GARCH 
or SV and their variations, which have found applications in many different areas. In general, 
the optimal predictor minimizes the expected loss function that measures a distance between the 
predicted and the actual value of the series. The most popular is the L 2 loss function since it is 
easily tractable and intuitively clear. In that case, the best predictor, in the sense of minimizing 
the mean squared error, is the conditional expectation. The mean squared error indicates the 
uncertainty in the prediction. However, this measure is not applicable when concerning random 
elements with no finite second moments, which are often observed in the reality. 

The same problem arises in the parameter estimation where the quality of estimators is tradi- 
tionally measured by the mean squared error. The estimator directly reflected such a criterion is 
the ordinary least square estimator (OLS) in linear regression analysis, which minimizes the L 2 
vertical distance between the observed data and the assumed linear model. A large number of 
related estimators have been invented and compared with OLS by means of L 2 distances between 
the true parameter and proposed estimators. However, once we lose the finite second moment 
condition for random components of the estimators, we could not use the mean squared error to 
judge the goodness of estimators. Then alternative measures are required. 

Although there are many potential loss functions, the LP loss function with < p < 2 is a natural 
plausible candidate to evaluate the goodness of prediction or estimators. Thus, we have to study 
the fractional absolute moments m p := E[\X\ P ] of order < p < 2. We also consider /i-centered 
moments mn,p '■= E[\X — fx\ p ] with fx € R. As far as we know, there have not been enough researches 
of the fractional absolute moments except for the special case of first order absolute moment m\. 
The reason is that the existing methods are unfamiliar or these methods seem to require a lot of 
numerical work. 

Taking this into consideration, we summarize existing methods for obtaining the fractional abso- 
lute moments. In particular, we focus on the methods exploiting the Laplace (LP) transform or the 
characteristic function (ch.f.) of the corresponding distribution. It is well-known that the moments 
of integer orders are related to the derivatives of ch.f. or LP transform at zero. More generally, 
the theory of fractional calculus can be utilized in order to obtain the non-integer real moments. 
There are several works giving the relation between the fractional moments and the corresponding 
ch.f or LP transform. We refer to von Bahr [23], Brown [2], Ramachandran [18], Kawata [8j Sec. 
11.4], Wolfe [251 [261 [271 [28], Laue [El EE], Zolotarev [291 Sec. 2.1] and Paolella pH Sec. 8.3]. The 
methods using moment generating functions have also been studied, e.g. by Cressie et al. [I] and 
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Cressie and Borkent [3]. Prom these references, we shall explain [8], [10], [26] and [3] which are 
convenient for our purpose. 

Our main tool is the fractional derivative of order 7 = k + A with /c G N, < A < 1, see e.g. \W\ 
Eq. (2.1)], 

■t /(*)(<)_ /(*)(„) 
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■ du, i G 



where / is a complex-valued function, /( fc ) is its kth derivative and T is the Gamma function. We 
are mostly interested in the fractional absolute moments m\ + \ = £'[|X| 1+A ] with < A < 1. For 
this reason, we will need the fractional derivative of order 1 + A at zero, 
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The construction of our paper is as follows. In the remainder of Section [IJ we make a brief survey 
on the relation between the fractional absolute moments and fractional derivatives using references 
cited above. Several convenient formulae are also derived. In Section [2] we apply the mentioned 
methods to the infinitely divisible distributions and examine their fractional absolute moments. 
Heavy tailed distributions, such as stable, Pareto, geometric stable and Linnik distributions, are 
considered. Especially, in Section [3] we pay attention to compound Poisson distribution that is 
popular in applications. In the final section, several applications are presented. The fractional 
errors of predictions concerning stable distributions are explicitly calculated. In addition, the 
estimation errors in regression models are evaluated by the fractional absolute moments in heavy 
tailed cases. 

1.1. Fractional derivatives of Laplace transforms. Let F be a distribution function (d.f.) of 
a non-negative random variable X. Its LP transform is defined as 



4>(t) :- 



-tx 
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t > 0. 







We characterize the relation between moments of X and the fractional derivative of 4>. Since (f> is 
only defined on the positive real axis, we define its fractional derivative of order l + A, < A < 1, 

as 

A <j>'(u) - <t>'(t) 
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We start with the following result, which is a modification of Theorem 1 in 
Lemma 1.1. Let < A < 1 and let (ft be the LP transform of the d.f. F{x) such that F(x) = for 



x < 0. Then mi + \ exists if and only if cj)'(0+) exists and the fractional derivative 
exists, in which case 
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Proof. Suppose that mi + \ exists, then (f)'(u) exists for all u > and is equal to — / °° xe xu dF(x) 
and 0'(O+) exists and is equal to — J^°xdF(x). We substitute this into the definition (|1.2|) of 



fractional derivative and use Fubini's theorem to see that 
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Conversely, the existence of 4>'(0+) implies 

roo 

-<f>\0+) > / xe~ ux dF(x) = -cj)'(u) 



for any u > 0. Hence, the reverse argument yields mi+x < °°- D 

Remark 1.2. (1) We applied the fractional derivative of order A to the first derivative <ft' , whereas 
in |261 Theorem 1] the \th fractional derivative of (p was differentiated. The latter approach leads 
to 

d r A r°° <t>(u) - </>(t) 
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which can be also useful representation of mi+\. The right-hand side can be written as 
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Since Lebesgue 's dominated convergence theorem enables to make differentiation under the integral 
sign, we recover the formula in Lemma 

(2) In Proposition 5 of Cressie and Borkent [3] , they apply the fractional integral to the derivative 
of moment generating function. 

1.2. Fractional derivatives of characteristic functions. We denote the ch.f. of a random 
variable X with d.f. F by 

/oo 
e itx dF(x), t€M, 
-oo 

and denote that for X — fi with fj, £ R by 

tp^t) := e- W/ V(i), t e M. 

As stated before, there are series of researches regarding the relation between the fractional de- 
rivative of ip and the fractional absolute moment. Among them we mainly work on the result by 
|10j referring to others [26], [8] as supplementary results. Since the original result by |10j is quite 
general and theoretical, we need to fill the gap between the theory and applications. 

Lemma 1.3. Let < A < 1 and let <p be the ch.f. of an arbitrary d.f. F. Then m\+\ exists if and 
only if 

1. Re / — 1 , . du exists, and 

Jo u + 

1 - Re<p(t) 

2. lim — -r exists. 

t^o+ t 1+x 

In such a case 

(1-3) m 1+x = — K- r-r X Re f * - du 

sm^r(l-A) J u 

and the fractional absolute moment with center \i (rrini+x = E[\X — ]) is given by 

1 A r 1°° e^X-u) , „ f°° e^ u f'(-u) , •> 
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Proof. We begin by applying equation (2.6) in |10] Theorem 2.2] with n = and h = tp^, which 
gives 



(1.5) 
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Since the first derivative of the ch.f. ip^ is 

ip'^u) = -ifxe- iu ^tp{u) + e _lu V (tt), u e 
the fractional derivative of order 1 + A at zero, by (jl.ip . becomes 
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and hence, by noting that Re<//(0) = 0, 
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Inserting this to (jl.5|l , we obtain the desired result (jl.4p . We may decompose the integral into two 
parts as in (|1.4j) since the existence of the integral 
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follows from | Rey?( — u)| < 1 and 
| Im ip{—u)\ < 



sm(—ux) \ dF(x) A 1 < / |ux| dF(x) A 1. 
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Taking // = 0, we get (jl.3p . The equivalence of the conditions 1 and 2 and the existence of mi + \ 
follows from part (b) of [10], Theorem 2.2]. □ 

Another expression of the fractional absolute moment of order < 7 < 2 is given in Theorem 
11.4.3 of [8], which is stated as the following lemma. 

Lemma 1.4. A necessary and sufficient condition for the existence of m~, < 7 < 2, is that 
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Remark 1.5. Equation ()1.7|) can also be found as (2.1.9) in [29] or (8.30) in |17| . in both cases 
with differently written constant in front of the integral and with a typo contained. We briefly see 
the equivalence of Lemma and Lemma \1.4\ A simple calculation yields 
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where we use Fubini's theorem in the last step. Hence, we obtain the first part of Lemma \1.4 
Moreover, due to the integration by parts, for x £ R we have 
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If condition (jl.6p holds, then by the Lebesgue dominated convergence theorem, 
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and by Fubini's theorem, 



1 f°° f°° xsinux , iri/ . 1 rRe^f-M) , 

■du<LF(a;) = - - / ; dn<oo. 



Thus, conditions 1 and 2 in Lemma 11.31 are satisfied. We can prove the converse in a similar 
manner. Finally, if m\ + \ < oo then 

1 — Re ip(u) 
lim t— t— — - = 

u->o+ u 1+A 

follows from (II. 6p . TTus, together with (jl.8jl and (|1.9j) : yields the expression (II. 3p in Lemma \1.3[ 
Moreover, Kawata [U Theorem 11.4.4] has obtained expressions for m 7 , 2 < 7, in the form of 

/"oo ^ 2fc 

- Q / -- {1+A) [l - Re p(t) + £ * c,( 2fc ) (0) 



dn, 



where £ € N is such that 2^ < 7 < 2^ + 2 and Cg is a positive constant depending on I. In 
other context, Wolfe [26j has derived techniques to obtain moments of any rational order from the 
ch.f. which are in different manner. It is well known, see e.g. [15], that the absolute nth moment is 
obtained from the identity 

1 f a sinxt 

— lim / dt = sign(x), x E R, 

7T a->oo J_ a t 

which yields, with applications of Fubini's theorem, that 

/oo />oo 1 />oo J 

\x\ n dF(x) = / x n sign(x)dF(x) = — — j / («) + ^ (-«)] — . 

-00 J— 00 ■^ 7ri J— 00 u 

Wolfe [26] generalized this relation for calculating m 7 , 7 G R. 

2. INFINITELY DIVISIBLE DISTRIBUTIONS 

In this section, we examine the class of infinitely divisible (ID for short) distributions, whose 
general definitions and many distributional properties are given by their ch.f. Many well-known 
distributions belong to this class and there are magnitude of applications in different areas (finance, 
insurance, physics, astronomy etc.). Here we work on the distribution without Gaussian part, its 
ch.f. is 

(2.1) tp(t) =exp\i5u+ / (e ite - 1 - itxl {a .<i } )z/(dx)), t € R, 

where 5 £ R is a centering constant and v is the Levy measure satisfying ^({0}) = and J" R (|x| 2 A 
1) v(dx) < 00. For more details on the definition and properties, we refer to Sato |21j . 

Although we can not calculate m 7 from density functions, because they are not available for 
most ID distributions, we can directly apply the fractional derivative to ch.f. and obtain fractional 
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moments. An advantage is that we can check the existence of fractional moments by the Levy 
measure of ID distributions and we need to check neither condition 1 nor condition 2 of Lemma 
11.31 The following result is well-known criteria for moments (see e.g. [214 Corollary 25.8] or \2A\ 
Theorem 2]). In our case of interest -E[|X| 7 ], < 7 < 2, we find a simple proof and give it in 
Appendix. 

Lemma 2.1. Let X be an ID distribution with Levy measure v. Then for < 7 < 2, m 7 < 00 if 
and only if 

/ |x| 7 v(dx) < 00. 

J\x\>l 

In what follows, we present the examples. 

2.1. Stable distributions. As a representative of heavy tailed distributions we firstly consider 
stable distributions. A random variable X has a stable distribution with parameters < a < 2, 
0- > 0, -1 < p < 1 and 5 G R if its ch.f. has the form, cf. [2Q1 Definition 1.1.6], 

(2.2) <p(t) = exp{iSt - a a \t\ a oo{t)}, t € R, 

where 

'\ -i/3 tan 3p sign (t), if a ^ 1 

1 + i/3| sign(i) log \t\, if a = 1. 

It is well-known that if 7 < a < 2, the moment of order 7 exists, otherwise it does not exist, see 
e.g. [181 Sec. 4] or [20\ Property 1.2.16]. We briefly review the existing results on the moments. If 
< a < 1 and X is a stable subordinator with the LP transform given by E[e~ uX ] = exp{— a a u a }, 
then for —00 < 7 < a, 

which is shown by Wolfe |26l Sec. 4] or Shanbhag and Sreehari [22]. In symmetric case (f3 = 0) 
with 5 = 0, it is shown in Shanbhag and Sreehari |22l Theorem 3] that 

!0 A , 2 7 r((l + 7 )/2)r(l-7/q) 

(2 - 4) m "= r(l- 7 /2)r(l/2) " ' " 1 <T r <«. 

where they rely on the decomposition of the symmetric stable distribution (see also Section 25 in 
Sato |2T]). For general (5 and 5 = 0, the following relation is proved by two different methods in 
Section 8.3 of [Hj, see also [201 P- 18], 

where = (3 tan ^ and 



(2.3) w(t) 



(2.5) ?n 7 = K _1 r(l - -)(1 + 6 2 )^ cos -arctan 6< a 7 , -1 < 7 < a, 

a Va / 



_ jr(l- 7 )cos^, if 7^1, 

\f, if 7 = 1. 

Using the fractional derivative, we obtain from Lemma 11.31 not only another proof of (|2.5p but also 
formulae for fractional absolute /x-centered moments which seem to be new. 

Proposition 2.2. Let X have a stable distribution with real parameters a > 1, < 1, S = and 

a > 0. Then, for < A < a — 1, we have 



AT(1 - 
sin(4p)r(l - A)' 



(2.6) m 1+A = . TT ^ 1+A (1 + 



cos I (l — — ) arctan #| + 9 sin 

{(1-— ) axctanfl] 
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and for 

(2.7) m ul+x = r \fi / u-^+^e-^" sin Uu- 6a a u a ) du 

V ; M+ sin(^)r(l-A)l F 7o 1 ' 

+ aa a u a - x - 2 e- aaua cos (fiu - 8a a u a ) - 6 sin (fiu - 6a a u a ) d-u} , 
where 9 = /3tan If X is symmetric (/? = 0) ; i£ follows that 



Ar(i-^) 

sin(^)r(l-A)' 



(2-8) ^i+A = .Ta^w-," 

and 
(2-9) 

m M+A = rn-( 1+A )e-"%in(^)dn + a f°° u^' 2 ^ cos (^) d u\. 

M, + sin(4f )r(l - A) \a J K a > J V a > J 

Proof. We begin with the expression of m^i+A- Let <^(«) be the ch.f. of a stable distribution with 
5 = and a > 1. Since we have, for u > 0, 

Ime i/M V(-u) = exp(-a a u a ) sin (fiu - 9a a u a ) , 
Ree^V (--") = aa a u a ~ l e^{-a a u a ) cos (fiu - 6a a u a ) 

- a6a a u a - 1 exp(-a a w a ) sin (fiu - 9a a u a ) , 

inserting these into (jl.4p of Lemma 11.31 we get (|2.7p . For rnj^, we let /i = in (j2.7j) and use 
change of variables theorem to obtain 

A(J 1+A /-oo /-oo 

mi_i_A = t i / it <* e cos6>udu + # / it « e sin#udu>. 

1+A sin(f )r(l-A)Uo Jo J 

Now applying the formulae 

/ x a_1 e _6:c sin cx dx = — -= — ^-2: — 777 sin { a arctan - 1, [Rea > — 1, Reo > I Imcll, 
Jo (b 2 + c 2 ) a / 2 I b) 

f°° L(a) r c ~\ 

/ x a_1 e _b;r cos cx dx = — , „ cos < a arctan -> , [Rea > 0, Re6 > I Imcl], 

Jo (y + cr) -' 2 I 0) 

see (3.944,5) and (3.944,6) in [6], we get (j2.6j) . Finally, letting /3 = and applying change of 
variables, the symmetric case is obtained. □ 

After some manipulation one can show that (|2.8p coincides with f|2.4|) and (|2.6p coincides with 
(|2.5p for 7 = 1 + A. Figure Q] shows the fractional absolute moments with center fi, computed 
numerically from the representation (|2.7p . We remark that even if this representation includes 
some integral expressions, it would be useful since most stable distributions have no explicit density 
functions. 

2.2. Pareto law. Another heavy tailed distribution is the Pareto distribution which has density 
and ch.f. given by 

f(x) = a(l + x)- a -\ x>0, 

"OO 



poo 

<p(t) = a ^{l + y)- a - l dy 7 t€l, 
J 
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Figure 1. The moments m^i+x of stable distribution with parameters a = 1.8, 
(3 £ [—1, 1], 5 = and a = 1. We choose A = 0.5 and depict the dependence on \i 
for three choices of j3 (left) and the dependence on j3 for three choices of \i (right). 

respectively, with real positive parameter a > 0. This distribution belongs to ID distribution (see 
Remark 8.12 in |21j). The fractional absolute moment mx+\ exists if and only if 1 + A < a. Though 
the density function is explicit, we obtain m^^ + \ from the fractional derivative of ch.f. Since 

/>oo 

Ime i/M V(-u) = a / sin{(/x - y)u}(l + y)~ a ~ l dy, 
Jo 

poo 

Ree i ^<p'(-u) = -a / y sin{(/x - y)u}{\ + y)""" 1 dy, 
Jo 

due to (|1.4p of Lemma [L3l we have, for 1 < 1 + A < a, 

m^ 1+x =«{(//+ l) l+x - a B{a - 1 - A, 2 + A) + |—^- 2 Fi(l,a + 1,3 + A; -/i)|, 

where we use formulae (3.196,1), (3.196,2) and (3.761,4) in [6] together with Fubini's theorem. 
Here, B is the beta function and 2F1 is the Gauss hyper geometric function. 

Although the following examples are not always in ID distributions, they are closely related and 
could be heavy tailed. 

2.3. Geometric stable law. A geometric stable distribution has similar properties to the stable 
distribution. The ch.f. is given as 

<p(t) = [1 + a a \t\ a uj(t) -iSt]- 1 , ten, 

where < a < 2 and uj(t) is defined by (|2.3p . However, its density function has no analytical 
expression. The tail behavior is the same as that of stable distribution, see e.g. Kozubowski et al. 

Lemma 2.3. Let X has a geometric stable distribution with 5 = 0. Then, for 1 < 1 + A < a and 
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a -X-2 {(1 + a a u a ) cos fiu + 6 a a u a sin liu} {1 + a a u a + 6 2 a a u a } , 
+ 2aa / u t——. — - , n , 77, - du 



{(1 + a a u a ) 2 + (9a a u a ) 2 } 2 
and 

Xa 1+X r°° _i_+x {l + v) 2 + {6v) 2 + 26 2 v 

mi+A " s in(^)r(l-A)7o 77 
w/iere 6> = /3 tan ^ . 

Proof. The proof is an application of Lemma 11.31 to the following formulae for u > 0, 

. . 1 + a a u a - Wa a u a 
<p{-u) 



{l + a a u a ) 2 + (9a a u a ) 2 
and 

„ , aa a u a ~ 1 2(1 + a a u a ){(l + a a u a )aa a u a - 1 + (9a a u a )a6a a u a ~ 1 } 

if ( — U) = —— 775 77 ^7 + 



+ i 



1 + a a u a ) 2 + {6<j a u a ) 2 {(1 + a a u a ) 2 + (<9<7 a u Q ) 2 } 2 



(1 + a a u a ) 2 + (0o- Q ti a ) 2 {(1 + a a u a ) 2 + (6a a u a ) 2 } 2 

□ 



If we put 9 = the results coincide with the standard Linnik law case. 

2.4. Linnik law. We consider a version of Linnik distribution given by Linnik p3]. Its density 
function is not explicit, while its ch.f. has the form 

(p{t) = (i + a a \t\ a )-P, ten, 

where < a < 2 is the stability parameter, a > is the scale parameter and /3 > 0. By the method 
of fractional derivative, we recover the result of Lin [13] as 

A/3cx 1+A uh 1 + A 1 + A, 
sm(^)r(l - A) a a 

where 1 < 1 + A < a. The fractional absolute moment of order 1 < 1 + A < a with center \i S M. is 

A(J 1+A ,00 „-(l+A) sin £^ roo n a-A-2 cog 

2.10) m ul+x = 7 {- — -^-du + a/3 / -g-f- du }. 

V ' M, + sin(^)r(l - A) U 7o (1 + w")/ 3 (l + n )^ 1 J 

For (3 = 1 these equations coincide with those in Lemma 12.31 for 9 = 0. 

2.5. Combination of stable law and Linnik law. Since 

lim (1 + a a \t\ a /[3)-P = e-^l*!" , t6R, < a < 2, 

/3->oo 

a symmetric stable distribution is a limit of Linnik-type distributions. We consider their combina- 
tion keeping both exponents a to be identical. Let X be a symmetric stable with ch.f. cp(t) = e - '*'" 
and let Y be a Linnik-type distribution with ch.f. ip(t) = (1 + lil"//?) - ^. Then we may express 
E[\X — Y\ 1+x ] by taking expectation of (|2.10p with \x replaced by X and a = j3~ l / a . As the result 
we obtain 

E\\X -Y\ 1+x ] 



A/3 1 " 1 ^ 



sin(4f)r(l-A) 

1+X 



UOO /* OO 



A/3 1 "— r(l-i±A) r 1 + A n 1 + A , , 1 + A i 1 + A 



U 1 



,2-/3-^)8) +17(1-1 + ^,1-/3-1^;^)}, 



sin(^)r(l - A) l v a 7 ' a" /v a a 
where ?7 is the confluent hypergeometric function (9.210,2) of [6]. 
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2.6. Subordinator. For practical purpose, it is desirable to express mi + \ through Levy measure 
v since ID distributions without Gaussian part are completely characterized by centering parameter 
S and Levy measure. However, in the light of (|2.ip . such expressions seem to be too formal and too 
complicated, thus they seem to be not very useful. Here, we confine our interest to some well-known 
distributions. However, for small classes of ID distributions general expressions of mi+A by v are 
worth considering. We pick out the class of subordinator (positive valued ID distributions) and 
that of compound Poisson distributions, the latter is treated in Section O 

For subordinator, we apply Lemma 11.11 and obtain a relatively simple expression. The LP 
transform of a subordinator can be found in Theorem 30.1]. 

Proposition 2.4. Let X be a positive valued ID random variable with shift parameter 5 > and 
Levy measure v such that Jj ^(1 A \s\) z^(ds) < oo. The LP transform is given by 

^(t)=e*(-* ) , t>0, 



*(t) = 6t+ ( (e st -l)v(ds). 

i(0,oo) 



A r r° 1 - e*(- u ) , f°° f°° 1 — e - us e*(-") , i 
Jo -u^ dU + J S ^ ds) l dn /- 



where 



Then it follows that 

3. Compound Poisson distribution 

Among ID distributions we focus on compound Poisson (CP for short) distribution which can 
easily manage the tail behavior by assuming a heavy tailed jump distribution. However, since most 
distributions do not have explicit representations, we rely on the ch.f. or the LP transform for 
calculating fractional moments. Let c be the intensity parameter of underlying Poisson distribution 
and v jump measure. The CP distribution has the following ch.f. 

(3.1) V9(t) = exp|c(y'(e to -l) l /(dx))| := exp{c(pj(t) - 1)} , t € R, 

where ipj(t) '■= J e ltx v(dx) is the ch.f. of jump distribution. If the jump distribution has positive 
support, we obtain the LP transform 

(j)(t) := exp{c(0j(t) - 1)}, *>0, 

where 4>j(t) := J e~ tx v(dx). The fractional absolute moments are expressed in the following lemma. 
The proof is just an application of Lemma 1 1.1 1 and Lemma 11.31 and we omit it. 

Lemma 3.1. Let tp(t) be the ch.f. of CP given by (|3.ip . then we have the following form for 
fractional fi-centered moments of order 1 < 1 + A < 2, 

A r f°° 

\i \ u - ( 1+A ) sin{/iu + clm(ipj(— u))} exp{c(Re(v9j(— u)) — 1)} du 
L Jo 



m 



sin(^f )r(l - A) 

/■OO 

+ c / u~^ l+x ^ Re((p'j(—u)) cos{fj,u + clm(<^ j(—u))} exp{c(Re((/;j(— u)) — 1)} du 



oo 

(1+A) 



u ( + >lm((p'j(—u))sm{fiu + cIm.(ipj(—u))}exp{c(R.e(<pj(—u)) — l)}du 

If the jump distribution is symmetric, i.e. Im <pj(u) = 0, we have 

A 



sin(f )r(l - A) 



fi u ( 1+A ) sm(fj,u) exp{c(tpj(— u) — 1)} du 
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/•OO 

+ c / u - ( 1+A ) ip'j(-u) cos(fiu) exjp{c(ipj(—u) — 1)} du 
Jo 

and moreover 



\r f'°° 

Xn , rn rr / u-^^'A-u) exp{ C (^(-u) - 1)} du. 
-75-11 (1 - A) Jo 



mi+X sm(f )r(l - A) 

If the jump distribution has positive support, we have 

\r f'°° 

mi+x = fxT^A) I u - (1+x) W u )^p{c(Mu) - i)} -&(o)}d«. 

In what follows, we will examine jumps given by well known distributions, which are not always 
heavy tailed, and try to obtain analytical expressions. Since they require a lot of numerical integrals 
and special functions we do not write complete calculation processes and just mention the key steps 
of derivation. 

3.1. exponential jump. The LP transform of the exponential distribution with parameter (3, 
i.e. with density function ^e x ^, x > 0, is <j)j(u) = 1/(1 + (3u), u > 0. Then due to Lemma 13.11 
fractional absolute moments for < A < 1 are given by 

Ac/3 f°° (1 + M 2 ~ exp{c(ii " 1)} 



mi+A = FU^A)/ + dU 

= c/3 1+A r(2 + A){iFi(l - A; 2; -c) + ^F X {1 - A; 3; -c)}, 

where \F\ is the confluent hypergeometric function (9.210,1) of [6] and we use (3.383,1) and (3.191,3) 

of m- 

3.2. symmetric stable jump. Recall that the ch.f. is (j)j(u) = e~' u ' Q with 1 < a < 2 and thus we 
apply Lemma 13. 1 1 with l<l + A<ato obtain the following series representation, 



\ac 


sin(^)r(l 


-A) 


Ac 




sin(4f)r(l 


-A) 


Ac 





mi+\ = — 7^ / u a - x - 2 e~ ua exp{c(e- uQ - 1)} du 

N Jo 

/•oo 

/ v - {1+x)/a e- v exp{c(e- v - 1)} dv 
Jo 

c yVn r(l-(l + A)/q) 
sin(^)r(l-A) ^ nl^ + l) 1 "^)/"' 

Again by Lemma 13. 11 shifted fractional moments E[\X — fi\ 1+x ] with 1 < 1 + A < a are obtained as 



A 



/•OO 

y. \ u~ {1+x) sin(/m) exp{c(e-" Q - 1)} du 
Jo 



sin(f )r(l - A) 

/•oo 

+ ca u a ~ x ~ 2 cos(^u)e- ua exp{c(e-" Q - 1)} du 
Jo 

3.3. Linnik distribution jump. Let ip j be the ch.f. of Linnik distribution with parameters a > 1, 
f3 > and a = 1. Since 

ip'j(-u) = a/3(l + li")-^ 1 /" 1 , u > 0, 
from Lemma 13.11 and change of variables formula (v = (1 + u a )~^) it follows that 

Ace~ c f 1 i±* i x _i±A „, , 

mi + \ = r / v <*p (1 — vf>) « e dv 

1+A srn(^)r(l - A) Jo V ' 
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for 1 < 1 + A < a. If (3 = 1, the jump distribution is the symmetric geometric stable distribution 
and we have 

Ace- C 1 + A 1 + A 1 + A 

m i+A = . tx^t: — rr^t 1 > x + Ji^n 1 + ; 2 ; c )> 



sin(^)r(l-A) a 



a a 



where we use (3.383,1) in [6]. 

3.4. deterministic jump of size 1 (simple Poisson). Substituting its LP transform 4>j(u) 
e~ u into the expression in Lemma 13-H we have 

Ace- C f°° e c - e~ u e cc ~ u , 
m i+A = FTi TT / ~TTx du > 



r(i - A) J 

which is rewritten by the Taylor expansion as 



X ce -c ™ck f oo 1 _ e -(k+i)u « k i+x 



c 



T(l - A) ^ k\ Jo u l + x ^ kl ' 

v ; k=0 JU fc=0 
where the final expression can be directly obtained from probability mass function. 

Remark 3.2. If the jump distribution has reproductive property, i.e. it is convolution- closed, we 
have another method for determining the fractional absolute moments. Write the CP random 
variable as SV = ^2jL\Xj, where N has the Poisson distribution with parameter c and (Xj) is an 
iid sequence such that X\ has reproduction property. Denote the ch.f. of kth convolution of X\ by 
<Pk(u)> then under suitable conditions we have 

A 



toi+a = t E Re / — du 

+ sin(^)r(l - A) L 



A 



00 Jt 



sin(f )r(l - A) § 



k\ 





1 P , n(- u ) 


Jo 


u i+x 




poo 


e" c 


Re / 




Jo 



fk(~ u ) 

u i+x 



du 



In case of the LP transform we denote that of kth convolution of X\ by <pk( u ), u > 0, and from 
Lemma \l.l\ we obtain 

A r r ^-j^(0+) I Ae- ^ c k ^(u) - ^(0+) 

4. Applications 

4.1. Evaluation of conditional expectation for stable law. Conditional expectations of stable 
random vectors have been intensively investigated in [7] and [19], since stable laws are often thought 
as natural generalization of the Gaussian random vector for which the minimizer of the mean 
squared error given some components of the vector is the conditional expectation. However, their 
evaluations have not been examined enough. In what follows, we evaluate the goodness of several 
predictors given by conditional expectations through their fractional moments. 
Firstly, we consider general results for a bivariate stable random vector with ch.f. 

y{t x M) ■■= E[e^ XM % (h,t 2 )eR 2 , 

which can be written as 

(4.1) (p(h,t 2 ) = exp { - J \t lSl + t 2 s 2 \ a (l - itan ™ signal + t 2 s 2 f) L(ds) + ifa^ + t 2 5 2 )}, 

where T is a finite measure on the unit sphere S 1 , called spectral measure, and we let a > 1, see 
[20\ Theorem 2.3.1]. Our aim is to linearly approximate X 2 by X\ and evaluate the fractional 
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error of order 1 < 7 < 2. The situation includes various settings, e.g. if stable random vectors are 
symmetric, i.e. 



<p(h,t 2 ) = exp I - J ^ \t\s\ + t 2 s 2 | Q r(ds)}, 



then it is proved that E[X 2 \ X{\ = cX\ with some constant c, see |20| Theorem 4.1.2] or |19t 
Theorem 3.1]. For general case we refer to [7J Theorem 3.1]. For convenience, we assume (6\, 6 2 ) = 
0, the general result for (6±,6 2 ) 7^ can be obtained in the same manner. 

Proposition 4.1. Let (X\,X 2 ) be a bivariate stable random vector defined by (|4.ip such that 
(61,62) = 0. Then for any constant c and l<l + A<a, it follows that 



Ar(i-i±A) 



E[\X 2 -cX x \^\ 



sin(f )r(l - A) 



l+A 1 



cos I (l — — ) arctan #o| + #0 sm | (l ~~ — ) arctan #o| 



1 + A N 



l/a f sl sign(s 2 -cs 1 )\s 2 -cs 1 \ a T(ds) 
s 2 -csi\ T(ds)\ , (3 = — 



where 

(4.2) ao = 

and 9q = /3q tan ^ . In symmetric case, we have 
E[\X 2 -cX 1 \ 1+x ] 



J s i \s2 ~ csi| Q r(ds) 



Ar(i 



l+A - 

a ' 



sin(^)r(l-A) V7 s i 



\s 2 — cs\\ a V(ds 



l+A 



Proof. The fractional derivative of the ch.f. of X 2 — cX\ is calculated. We put t\ = —cu and t 2 = u 
in (|4.ip . then we regard it as a function of u, 

^! e iu{X 2 -cX 1 h 



exp 



{-\< [ 



\c*r(A \fi ■+ 7ra • ( x/§isign(s 2 -csi)|s 2 -csi| a r(ds) 

s 2 -cs x \ r(ds) 1 -itan — sign(n) =—j — — 

V 2 J gl |s 2 - csi\ a T{ds) 



)}• 



In view of (|2.2p . this is the ch.f. of one-dimensional a-stable distribution with parameters (/3, cr, 5) = 
(/?o,a"o,0). Hence, we apply Proposition 12.21 to obtain the result. □ 

Examples. As examples we consider predictions for two bivariate stable random vectors and one 
stable process. First we treat a bivariate stable random vector considered by [TBI p. 183] such that 
ch.f. of (X\,X 2 ) satisfies for \a\ < 1, a / 1, 



¥>(ti, t 2 ) = E[S l ^ +t ^] = exp{ - a a \t 2 \ a (l -f i/3tan ^ sign(i 2 ) 

— sign(ti + at 2 ) 



1 



— I — j — 1 1\ + at 2 \ a [l + i/3 tan 



2 1 — sign(a)|a| 



The conditional ch.f. is 



<fXi=x(u) ■= E[e luX2 \X 1 = x}= expjkmt - a a \u\ a (l + i/3 tan — sig 

and hence for 1 < a < 2, £?[AT 2 | Ali = x] = ax. The support of spectral measure T consists of four 
points in S , 



Tra 



n(u))} 



r(0,±l) = -^(l±/3), 



r ± 



vTT^' VT+~o?J 21 -\a 



1 <7 C 



: (l + a 2 )f l±/3 



1 — sign(a)|o| 
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Hence, 

la — c| 



f |s 2 -csi| Q r(ds) = a a (l + 

/ sign(s 2 - c Sl )\s 2 - c Sl \ a T(ds) = pa a (l + 
Js 1 v 



1 - \a\ a . 
sign(a — c)\a — c\ 



Substitution of these relations into (|4.2|) yields o"o and /3q that can be used for calculating the 
fractional absolute prediction error by Proposition 14. 11 

Another example is the prediction for sub-Gaussian random vector. Let < a < 2, |7| < 1, and 
let (G±, G2) be zero mean Gaussian random vector with covariance matrix 

'1 ^ 
,7 1 



(4.3) 



Let A be a positive a/2-stable random variable, given by the LP transform 

E[e~ uA ] =e" W2 , u>0, 

such that it is independent of (Gi,^). The vector {X\,X2) = ( y A 1 l 2 G\, A X I 2 G2) is called a sub- 
Gaussian symmetric a-stable random vector. In [19j, E\X^ \ Xi] = jXi is shown. Since we have 
the ch.f. 

^ [e m(X 2 - 7 Xx)] = E[E[e iu(A^G 2 -yA^G 1 ) ] | ^ = ^ [exp{ _^(_ 7j l) S (- 7 , 1)'}] 

= £[exp{ { - Y - LlA ^ = eM-(^ L T /2 u a }, 

due to the fractional moment (|2.8p . we get the fractional error 



E[\X 2 -E[X 2 I X 



sin(^)r(l - A) V 2 



where 1 < 1 + A < a. If A 2 is predicted by a linear function cX±, in a similar manner, we obtain 

„2 \ a/2 



E[e ^-^)] = exp - ( 1 " 27 2 C + C ) „« 



which yields 

11 J sin(f )r(l - A) V 2 ; 

Alternatively, we could use the spectral measure of sub-Gaussian random vector given in \20\ 
Proposition 2.5.8]. Since it is given in a closed form, we obtain the fractional error directly from 
ch.f. here. 

Next we examine the prediction of the a-stable OU process with < a < 2 and 7 > given by 

X t = e-^Xo + f e"^- s ) dZ s , t> 0, 
Jo 

where {Zt}teR is the symmetric a-stable motion. We set Xq = e 7S dZ s to obtain the stationary 
version, see [20|, Example 3.6.3] for its definition. Then, the conditional ch.f. of Xt given Xq is 

ip Xo (u) = E[e luXt I X ] = exp {me-^Xo} exp |- J |ue" 7 (*- s )| a da j 

f 1 - e~ Q7 ' ] 
= exp < iue 7 X |ti| a > , 

I «7 J 
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which yields E[X t \ Xq] = e _7 *Xo for a > 1. Since the mean squared error of the prediction is 
not available, we use the fractional absolute moment of order 1 < 1 + A < a. More generally, we 
measure the error of a linear approximation cXq with c G M. 

Proposition 4.2. Let Xt be an a-stable OU-process driven by the symmetric stable motion with 
the location parameter 5 = 0. Then for 1 < 1 + A < a, 



EUXt-cX, 



o 



'(1-^) f 
sin(4f)r(l-A)l 



Ar (l_i±A) f J _ e -ajt + ( c _ e - 7 t)a . 1±A 



«7 



and hence putting c = e 7 , we obtain 



E[\X t -E[X t | X 



J sin(4nr(l - A) V «7 ' 



Proof. Since 



E r e iu(X t -cX ) | Xq \ = eX p{ ilt ( e -7t _ c )x }exp{- — (1 - e- Q7i )|n| a }, u G 

we may use formula (|2.9|) in Proposition 12.21 with /i = (c — e _7 *)Xo and (j™ = ^(1 — e 
Consequently, 



E\\X t - cX \ 1+x I X 0J - 



A(J 1+A r (c - e" 7 *)^ 



r(i-A)l 



sin(4f )r(l - A) 



sin 



+ a 



(T 



v a ~ x ~ 2 COS 



(c - e" 7 ')^ 



(c - e" 7 *)X 



(T 



e""" dv}. 



- v " dv 



After taking expectation w.r.t. Xq and applying Fubini's theorem, we use 



Ecos 



(c - e~^)X 



a 



E 



(c-e-^Xo . i(c-e-^)X 



sm 



a 



a 



exp 



1 /c-e" 7 * 
«7 



c — e 



(7 



a 



1 /c-e" 7 * 



exp <^ 

07 



(7 



to get 



£[|JQ-cXo| 1+A ] 



l+A 



sin(^)r(l - A) 



1 + 



1 / c — e 7 *\ a 



07 V O" 



l+A 

it « exp 



07 \ <7 



M > dit, 



where we apply change of variables several times. Then the result is implied by a a = — (1 — e Q7 *) 
and definition of the gamma function. □ 



Since the finite dimensional distribution of the a-stable OU process is a multivariate stable, we 
may use Proposition 14.11 similarly as before. The spectral measure is given in |19l Example 3.6.4]. 



4.2. Evaluation of conditional expectation for Linnik law. Let (Xi, X2) be a bivariate Linnik 
distribution with ch.f. 

(p^, ta) = {1 + (i'St) Q / 2 }-^, t' = (h,t 2 ) G M 2 , 
where < a < 2, /3 > and £ is given by (|4.3p . see e.g. [12] , 
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Proposition 4.3. Let (Xi,X2) be a bivariate Linnik random vector. Then E[X<2 \ X\\ = 7X1 and 

for any c S R and 1 < a < 2 it follows that 

E[\X 2 - cX^) = y U> - 2 7C + ll^fld - 1±*./, + 1±*) 

sm(— )r(l - A) 

and therefore 



a a 



2 

Proof. To obtain £7[X2 | Xi], we use the decomposition by [5] of univariate Linnik law, which is 
also applicable in our bivariate case. Let (Yi, Y 2 ) be a sub-Gaussian random vector with ch.f. 

v(t l ,t 2 ) = e-( t '^ a/2 , t = ( tl ,t 2 y 

and let Z be an independent random variable with density 

f{x) = wwr x>0 - 

Then we observe that (X 1 ,X 2 ) = (Y 1 Z 1 / a > 3 , Y 2 Z 1 l a P), which leads to 
E[X 2 I Xi] = ElE^Z 1 ^ I y^Z 1 /^] I YiZ 1 ^} 

= E[z x i afi E^ 2 1 y] I y x z x i^\ = 7 yiZ 1 /°' 3 = 7X1, 

where the conditional expectation of the sub-Gaussian random vector is used. Now put t\ = —cu 
and t 2 = u in tp(ti,t 2 ) to obtain 

E ^iu(X 2 -cX 1 ^ = {1 + ( c 2 _ 2 7C + l)"/ 2 ^!"}-' 3 

and we conclude our result from Subsection 12.41 

□ 

4.3. Estimation errors of regression model. We consider the basic regression model 

Yj = 6 + Xj9t + £j, j = 1,2, . . . ,n, 

where (ej) is an iid sequence of symmetric random variables. It is well-known that the least squares 
estimator (9q, 9%), which is the best linear unbiased estimator if the e,- follow Gaussian distribution, 
has the form 

6 = Y - 9ix, 9 X 



E^ =1 (x J -x)(y J -y) 



Y%=i(xj - x) 2 

where x = - Y^j=\ x j an d Y = - Y^j=i Yj- For our purpose, it will be convenient to rewrite this as 

2- " fa - x)x - E]=i(xj ~x) 2 /n 

"0 — "0 — — ; =:r7; e,: ' 



8=1 





^E^iixj ~ x) 2 

We express the fractional errors for the case of a-stable noise distribution. In a similar manner, 
it would be possible to calculate the fractional errors for other regression- type estimators, e.g. [1], 
and compare the goodness of estimators. 

If £1 is a standard symmetric stable random variable with parameters 6 = 0, a = 1 and a > 1, 
then the characteristic functions of estimation errors are 

^«(S**-«*)] =e -^H" j = 0,1, 
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where 



2 



This together with f|2.8j) yields, for 1 < 1 + A < a, 



mOk-Bu\ 1+x \= [A -"", k = 0,l. 

sm(^-)r(l - A) 

Interestingly, if e = (e\, £2, ■ ■ • , £ n ) is an elliptically contoured stable random vector with ch.f. 
E[e [u ' E ] = e -\ u ' Iu \ a/2 for u' G R n , where I is n x n identity matrix, then we obtain closer results 
to Gaussian case. Namely, for 1 < 1 + A < a, 



E[\e k -e \^]= fc = o,i, 



l+Ai _ yyl V x a ) — 1+A 
sin(^)r(l - A) 

where 

x 2 1 \ 1/2 



°"0 = — 7 =7^ + ~ an d 0"! 



•>:; ,(•'•, *) 2 (E-=i(^-^) 2 ) l/2 ' 

Moreover, if £ = (ej, £2, . . . , £ n ) is a multivariate Linnik random vector with ch.f. ^[e 1 ' 
{1 + {u'Iu) a l 2 }-> 3 for u' G R n , we obtain in a similar manner that (1 < 1 + A < a) 



E[\e k -e k \ 1+x ] = xmi x 1 ^^ +1 ^ X) wi +x , k = 0,1. 

sin(^)r(l - A) 



Appendix A. Proof of Lemma 12.11 
Proof. We express ch.f. ip(t), given by (12.1j) . as the product ipi(t) • 922 (*) , where 

<pi(t) := exp iiSt + / (e ite - 1 - ltx) v(dx)\, t G R, 
ip 2 (t) := exp { [ (e itx - l)v(dx)\, t G R. 

W|a:|>l J 

Since the distribution with ch.f. ^i(i) has moments of any positive order, it suffices to consider 
<f2(t)- We use the necessary and sufficient condition (]1.6p for the existence of m 7 , see Lemma 11.41 
The following inequalities 

a 



C 1 - e~ a < a, a > 0, 



1 + a 

b 2 b 

1 - cos 6 < — < -, < b < 1, 



and the fact 

j (1 - cos fa) u(dx) < j ( ^ A l) u(dx) := c t < 00 

y|as|>l V|x|>l V 2 ' 

[ to obtain 

1 — Re ^2(i) > 1 — exp j A (cosfa — 1) z^(dx)| > J (1 — cosfa;) v(dx), 

L 7|x|>l J 1 + Ct J| X |>1 

1 — Re if2(t) = 1 — exp j J (cosfa; — 1) z/(dx)| + exp | J (cos tx — 1) z/(dx)| 

sin fa; i/(dx) ) exp < / (cos tx — 1) v(dx) > 
>i ' J|x-|>l - 1 
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< I (1 — cos tx) v(dx) -\ — / |sintx|^(dx) 

J\x\>i 2 J| x .|>i 

Then we notice by Fubini's theorem that 



t -(i+7) / (l-cos*x)i/(dx)dt= / rr dv / |xp i/(dx), 

J\x\>l J0 ^ +7 J|x|>l 



/ I sintxl u(dx) dt = / — r— — dt> / |x| 7 i/(dx). 

'|x|>l Jo « +7 ./|s|>l 

Since double integrals on the left-hand sides exist if and only if the integrals on the right-hand sides 
exist, the condition (jl.6p is equivalent to the existence of f\ x \ >1 |x| 7 ^(dx). □ 
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